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Abstract 

We study the spin-orbit splitting (Zs-splitting) for the proton d-orbits in 23 F in the Hartree- 
Fock framework with the tensor force in the exchange channel. 23 F has one more proton around 
the neutron-rich nucleus 22 0. A recent experiment indicates that the Zs-splitting for the proton d- 
orbits in 23 F is reduced from that in 17 F. Our calculation shows that the Zs-splitting in 23 F becomes 
smaller by about a few MeV due to the tensor force. This effect comes from the interaction between 
the valence proton and the occupied neutrons in the 0cZ 5 / 2 orbit through the tensor force and makes 
the Zs-splitting in 23 F close to the experimental data. 
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I. INTRODUCTION 



The spin-orbit splitting (/s-splitting) is important for the structure of nuclei. A large Is- 
splitting between single-particle orbits with the same orbital angular momentum is respon- 
sible for the shell structure of nuclei lj . Recently we have been obtaining much information 
about unstable nuclei from various experiments. There are experimental evidences, which 
indicate that the shell structure in neutron-rich nuclei changes from that in stable nuclei. 
To confirm the change of the shell structure, the information about single-particle orbits 
around closed-shell or sub-closed-shell nuclei is important. Michimasa and his collaborators 
studied the proton single-particle orbits in 23 F experimentally through the proton transfer 
reaction |2|. 23 F has one more proton around 22 0. They reported that the /s-splitting for 
the proton d-orbits_(5/2 + -3/2 + ) is 4.06MeV, while the /s-splitting for the proton (i-orbits in 
17 F is 5.00MeV which is similar to the neutron cZ-orbits in 17 (5.08MeV) 3, 4] due to 

the isospin symmetry. It indicates that there is a possibility that the /s-splitting is changed 
by the excess neutrons around 16 0. The shell model calculations reproduce the change of 
the /s-splitting from 17 F to 23 F nicely [2]. In the shell model calculation the /s-splitting in 
17 F ( 17 0) is an input parameter. Hence, it is interesting to study the /s-splitting with a 
mean-field-type model, where the /s-splitting is obtained self-consistently. 

Hartree-Fock and Hartree-Fock-Bogoliubov calculations can now be performed in the 
whole mass region over the nuclear chart. Such mean-field calculations can reproduce binding 
energies and radii of nuclei including unstable ones using effective forces with relatively 

nn 

simple forms like the Skyrme or Gogny force |5|, 16|]. In the mean field calculations the /s- 
splitting of single-particle orbits is produced mainly by the LS force. The /s-splitting of 
single-particle orbits and the magic number for binding energies can be explained with the 
LS force having the same strength in almost the whole mass region at least near the stability 
line. Some studies show that the /s-splitting in neutron-rich nuclei becomes small because 
;he diffuseness of the neutron density becomes large and the spin-orbit potential is weakened 

The tensor force acts on the spin of nucleon directly and should affect the /s-splitting. 
Although the tensor force is not usually included in the mean field calculations, some Hartree- 
Fock calculations explicitly including the tensor force or the pion in the relativistic model 
showed that the tensor force affects the /s-splitting in spin-unsaturated nuclei 
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Only one orbit of the spin-orbit partners is occupied in a spin-unsaturated 



nucleus, while both the spin-orbit partners are fully occupied in a spin-saturated nucleus. 
For example, 48 Ca is a spin-unsaturated nucleus, where the neutron O/7/2 orbit is a spin- 
unsaturated orbit and 40 Ca is a spin-saturated nucleus. Because the total spin coming from 
the intrinsic spin of nucleon is zero in a spin-saturated nucleus if the wave functions of 
spin-orbit partners have the same radial forms, the tensor force does not act between the 
spin-saturated core and a particle or a hole around the core. In a spin-unsaturated nucleus, 
the total intrinsic spin coming from the spin-unsaturated orbit has a finite value and the 
tensor force becomes active. In fact, the sizes of the /s-splitting for hole orbits change from 



40/ 



Ca to 48 Ca and lb O to 22 in t 



re results of the Hartree-Fock calculations with the tensor 



11 



12 



13 



14j. For the calcium isotopes, there is an experimental 



force or the pion jlO 

evidence 17| that the /s-splitting becomes smaller from 40 Ca to 48 Ca and the order of the 
change is comparable to that induced by the tensor force or the pion [lj[ Ijj]. It should be 
noted that in the Hartee-Fock approximation the energy contribution from the tensor force 
or the pion from the direct channel becomes zero and only that from the exchange channel 
has a finite value in closed-shell nuclei. 

Otsuka and his collaborators discussed the effect of the tensor force on single-particle 
energy in other mass regions. They nicely reproduced the change of the splitting between 
7r 0/iii/2 and TrOgj/2 in the Sb isotopes with neutron number 18| by the monopole shift 
induced by the tensor force 14j. They also suggested the effect of the tensor force on the 
shell evolution in the neutron-rich sd- and p /-shell region [l4,[l^]. They discussed that the 
neutron shell structure changes with proton number due to the monopole interaction between 
proton and neutron orbits and explained the appearance of the magic number 16 and the 
disappearance of the magic number 20 in the neutron-rich s<i-shell region [y , 19, 20]. They 
claimed that the monopole interaction is caused by the tensor force Q, lj|. To confirm 
such a discussion, the direct information about a single-particle state is essential. 

In this paper we perform the Hartree-Fock calculation for 22 and 23 F. We include the 
tensor force and study its effect on the /s-splitting. We also calculate 15 > 16 > 17 to see the effect 
of valence neutrons on the /s-splitting and its relation to the tensor force by comparing with 
22 and 23 F. The formulation is given in Section [Til and the results are given in Section [TTTI 
Section IIVI is devoted to the summary of the paper. 
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II. FORMULATION 



In the preset paper we adopt two types of Hamiltonian. One includes the 3-body force 
in addition to the kinetic term and the two-body force. The other includes the density- 
dependent force instead of the 3-body force. The Hamiltonian with the 3-body force H 3B 
and that with the density dependent force H BD have the following forms, 

H3B = J2^ + t, v ^ + E ^\n^r k )-E cm (i) 

i=l i<j i<j<k 

H ™ = E m + + E- (DD) (p;n,r,) - E CM . (2) 

i=l i<j i<j 

In the above expression, p, r, and M are the momentum, coordinate including spin and 
isospin, and mass of nucleon respectively. A is a mass number, v and are the 2-body and 
3-body potentials respectively. 

V (DD) ig 

the density-dependent potential with the one-body 
density p. We subtract the energy of the center of mass motion Eqm = P'i) 2 /2 AM. 

In the Hartree-Fock calculation we assume the wave function of the nucleus has the 
following form, 

V=Al[Mr a ) (3) 

a 

with the antisymmetrization operator A for nucleon coordinates, a labels each single-particle 
state and runs over all occupied state. With the wave function the total energies become 

E 3B = J2M^\ila) +Y,(^M V \^) + E (^a^> (3) |«^ 7 > ( 4 ) 
a ~ a<fS a</3<7 

for H 3B and 

£ DD = Z>«l|^«> + Z>«^MfiS> + J](^k (DD) (p)l^> (5) 

a ' a</3 a</3 

for i/ DD , where the tildes represent the antisymmetrization. In the above equations, E C m 
is dropped for simplicity. By taking a variation of the total energy with respect to a single- 
particle wave function ip a , we obtain the Hartree-Fock equation for each case: 
P 2 



2M 



1 Ha J J 1 

+ {^j(y)^a(z) - tp a (y)^(z)} ipp(x) + {tp a (y)^p(z) ~ ^p(y)^a(z)} ^ 7 (x) = e a ip a (x) (6) 



ip a (x) + E / d y^l(y) v ( x >y) [foivWaix) - ^Mfoix) 
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for the three-body force case and 
P 2 f 



p 



+ [ d y f dz ^l(y)ijl{z)—r — (p;y,z)TTf( x ) [M^Vvfc) -^M^p ( z )\ = e ^ a {x) 

~J J op 5^ a 



2M 

P ' 

6v^) f _ 6p_ 

(7) 

for the density-dependent force case. In the above expression the integrations over y and z 
include the summation over the spin and isospin index. 

In the present study we assume each single-particle state as an eigenfunction of total spin 
j = I + s. With the assumption a single-particle wave function can be expressed as 

Mr) = R«{r)y lajama (n)((jji a ), (8) 

where R is a radial wave function, y is an eigenfunction of j, and ( is an isospin wave 
function, a stands for node n a , total spin j a , its projection on the z axis m a , and isospin 
\i a . We do not assume the degeneracy for the orbits with the same n a , j a , and p a because 
the spherical symmetry of a mean field is broken in odd nuclei. It means that the states with 
the same n, j, and p but different m's are allowed to have different radial wave functions. In 
such a case we need to perform an angular momentum projection to obtain a wave function 
with a good angular momentum. The expectation value for the total angular momentum J 2 
with the wave function obtained in the Hartree-Fock calculation for a one-particle or one- 
hole state does not deviate from j u (j u + 1) largely (less than 1%), where j u is the total spin 
of the particle or hole orbit. It indicates the obtained wave function is almost an eigenstate 
of angular momentum. Hence, we do not perform the angular momentum projection. 
We approximate the density in a density-dependent force as 

a 

for calculational convenience. This expression is exact for a closed-shell nucleus with the 
spherical symmetry and should be a good approximation for a one-particle or one-hole 
nucleus with almost a spherical core. 

We exp and a radial wave function R a (r) by Gaussian functions with widths of a geometric 



series 



211 ] . We take 11 Gaussian functions with the minimum width 0.5fm and the maximum 



width 7fm for each single-particle state. The Hartree-Fock equation is solved by the gradient 
or damped-gradient method [22]. 



III. RESULT 

In this section we apply the Hartree-Fock method to 15 > 16 ' 17 > 22 and 23 F. We assume 16 
as a closed-shell nucleus up to the Op-shell and 22 O as a sub-closed-shell nucleus where the 
neutron Oc/5/2 orbit is fully occupied in addition to the occupied orbits in 16 0. For 22 there is 
the experimental evidence which suggests it has the sub- closed- shell structure of the neutron 
0^5/2 orbit 



231 ] . In the 15 O case, one neutron is subtracted from the neutron Opi/2 orbit or 
the neutron 0^3/2 orbit in 16 O. In the 17 O case we add one neutron in the 0d 5 / 2 orbit around 
16 0. We do not put a neutron in the (W3/2 orbit in 17 because there are no bound states in 
this configuration. In the 23 F case we add a proton in the 0d§/2 or 0d 3 / 2 orbit around 22 0. 

As for the effective interaction, we adopt the modified Volkov force No. 1 (MVl)[24j for 
the central part and the G3RS force 2lj] for the tensor part. We also include the Coulomb 
force. The G3RS force is determined to reproduce the nucleon-nucleon scattering data and, 
therefore, the tensor force in the G3RS force is the one in the free space. For the strength of 
the tensor force in the nuclear medium we do not have a definite guideline at present. The 
effective interaction obtained from the G-matrix th eory has a tensor part with a strength 
comparable to the tensor force in the free space Jj], 26, 27, 28] at least in the region where 
the relative distance is greater than about 0.8fm. We use the tensor force in the free space in 
the present calculation but we need a further investigation to determine the strength of the 
tensor force to be used in a mean field calculation. It should be noted that the difference in 



;he short range (r < 0.8fm) does not influence the tensor 



brce matrix elements significantly 



28| . As for the LS force we take the 5-type LS force 0, E]: 

iWofa + cr 2 ) • fc~ x <J(n - r 2 )~k . (10) 

The Majorana parameter in the MV1 force is fixed to 0.59, which is determined to reproduce 
the binding energy of 16 0. Wq in the LS force is taken as 115MeVfm 5 , which is the same as 
in the Gogny Dl force and is determined to reproduce the Zs-splitting for the Op orbits in 

In Table HJ the results for 16 0, 17 0, and 15 are summarized. The experimental data 



6 



TABLE I: Total energy (Etot), kinetic energy (T) and potential energy (V) of 16 0, 17 0, and 15 0. 
Vls and Vt are the contributions from the LS and tensor forces to the potential energy. Those are 
give in MeV. R c and R m are the charge and matter radii in fm. The last row shows the differences 
of energies between 15 (Op^) and 15 (Op^). In the parentheses the experimental data are 
given. 





E T Q T 




T 


V 


Vls 


Vt R c R m 


16Q 


-128.3 (- 


-127.6 a ) 


233.8 


-362.0 


-1.0 


0.0 2.71 (2.730(25) 6 ) 2.58 (2.54(02) c ) 


17 (0d 5/2 ) 


-132.3 (- 


-131.8 a ) 


254.7 


-387.0 


-4.1 


0.0 2.72 (2.662(26) 6 ) 2.64 (2.59(05) c ) 




-110.2 (- 


-112. a ) 


219.6 


-329.7 


-4.9 


-0.1 2.70 2.55 (2.44(04) c ) 


15 o (op-; 2 ) 


-104.5 




212.4 


-316.9 


0.9 


0.0 2.74 2.59 


A(0p 3 7 2 - 0p-/ 2 ) 


5.7 


(6.18 d ) 


-7.2 


12.8 


5.8 


0.1 



a Reference 2§| 
^Reference 3T| 
c Reference 
^Reference 



h 

32 1 



are also given in the parentheses if available. The potential energy from the tensor force 
becomes quite small because 16 is a LS-closed-shell nucleus. In the LS-closed- shell nucleus 
both the spin-orbit partners are completely occupied. Hence, the LS-closed-shell nucleus is 
a spin-saturated nucleus. The LS-closed-shell nucleus does not have a finite total orbital 
angular momentum and a finite total spin angular momentum. The tensor force consists of 
the rank 2 tensors of the orbital and spin angular momenta. Thus, the tensor force does 
not work between the LS-closed-shell nucleus and a particle or a hole around it, because a 
particle or hole has a spin angular momentum 1/2. In the last row the energy differences 
between 15 O (Op.^) and 15 O (Op^) are shown. It corresponds to the Zs-splitting for the 
Op orbits. It is about 10% smaller than the experimental value. The contribution from the 
LS force is 5.8MeV and is almost the same as the total Zs-splitting. It indicates that the 
/s-splitting is mainly produced by the LS force. The large contribution from the kinetic 
energy is almost canceled out with the contributions from the central and three-body forces. 
In 15 O the effect of the tensor force on the /s-splitting is negligible. 

In Table [TT] the results for 22 and 23 F are summarized. Although the binding energy 
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TABLE II: Total energy (-E7tot)> kinetic energy (T) and potential energy (V) of 22 and 23 F. Vls 
and Vt are the contributions from the LS and tensor forces to the potential energy. Those are 
give in MeV. R c and R m are the charge and matter radii in fm. The last row is the differences of 
energies between 23 F (0^3/2) and 23 F (0^5/2). In the parentheses the experimental data are given. 

^tot T V Vis Vt R c Rm 

22 -161.8 (-162. a ) 361.4 -523.2 -20.8 1.9 2.74 2.85 (2.88(06) b ) 

23 F (0d 3/2 ) -166.5 376.4 -542.8 -16.3 0.1 2.89 2.90 

23 F (0d 5/2 ) -170.7 (-175.3 a ) 383.9 -554.5 -24.1 3.2 2.84 2.87 (2.79(04) & ) 

A(0d 3/2 - 0d 5/2 ) 4.2 (4.06 c ) -7.5 11.7 7.8 -3.1 



a Reference 
b Reference 
c Reference 

of 23 F (Oc/5/2) (the ground state) is about 5MeV smaller than the experimental value, it 
probably does not affect our discussion on the /s-splitting. In 22 O the neutron 0^5/2 orbit 
around the 16 core is fully occupied. Because the spin-orbit partner, the neutron 0^3/2 orbit, 
is empty, 22 is a spin-unsaturated nucleus. Hence, 22 has a finite total orbital angular 
momentum and a finite total spin angular momentum, and the expectation value for the 
tensor potential energy in 23 F becomes finite. In 22 0, the energy contributions from the LS 
force and the tensor force are — 20.8MeV and 1.9MeV respectively. In 23 F a proton is added 
to 22 0. If the proton is put in the 0^3/2 orbit the absolute value of the LS potential energy 
becomes small by 4.5MeV and if the proton is put in the 0d§/2 orbit that of the LS potential 
energy becomes large by 3.3MeV. In contrast, the tensor potential energy becomes small by 
1.8MeV when the proton is in the 0<i3/2 orbit and becomes large by 1.3MeV when the proton 
is in the 0d 5 / 2 orbit. As a result, the contribution to the /s-splitting for the proton Od orbits 
in 23 F from the LS force is 7.8MeV and that from the tensor force is — 3.1MeV. The sum 
of them is 4.5MeV. The relatively small /s-splitting 4.2MeV after adding the contributions 
from the kinetic and other potential energies, which is close to the experimental value, is 
realized by the cancelation between the contributions from the LS force and the tensor force. 

The energy differences between one-particle states and their corresponding cores are 
shown in Table [TTTl The LS potential energies from the cores for the 0c?5/ 2 orbit are — 3.0MeV 
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TABLE III: Differences of the LS potential energy (A(Vls)) an d the tensor potential energy 
(A(Vt)) between one-particle nuclei and their core nuclei. Those are given in MeV. 



A(Fls) Agrj 

17 O(0d 5/2 ) - 16 O -3.0 0.0 

23 F(0d 3/2 ) - 22 O 4.5 -1.8 

23 F(0d 5/2 ) - 22 O -3.3 1.3 



in 17 and — 3.3MeV in 23 F. The LS potential energy for the Od^/2 orbit in 23 F is smaller 
as expected from that for the 0d 5 /2 orbit (3.3 x (2 + l)/2 « 5.0MeV). It is probably due 
to a weak binding of the Od^/2 orbit compared to the 0d 5 / 2 one. The contribution from the 
tensor force to the splitting for the Od orbits in 23 F is about a half of that from the LS force 
with the opposite sign as discussed in the previous section. The results for 17 and 23 F in 
Table III 1 1 indicate that the contribution to the /s-splitting from the LS force mainly comes 
from the 16 O core and that from the tensor force comes from the excess neutron orbit (the 
neutron 0d 5 / 2 orbit). 

TABLE IV: Potential energy contributions from the triplet-even tensor force (V^ E ) and the triplet- 
odd tensor force (VjP) in MeV. In the last two rows, the differences between 23 F (0^3/2 or 0^5/2) 
and 22 are given. 

y T 3E y T 3 ° 
22 o 0.1 1.8 

23 F (0d 3/2 ) -1.3 1.4 

23 F (0d 5/2 ) 1.0 2.1 

A( 23 F(0d 3/2 ) - 22 O) -1.4 -0.4 

A( 23 F(0d 5/2 ) - 22 O) 1.0 0.3 



In Table IIVt the contributions to the tensor potential energy from the triplet-even and 
triplet-odd parts are shown separately. In 22 O the tensor potential energy mainly comes 
from the triplet-odd part. It is natural because only the neutron 0^5/2 orbit is occupied 
and there are no valence protons around the 16 core. In 23 F the contribution from the 
triplet-even part is comparable to that from the triplet-odd part for the 0c?3/2 orbit and they 
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have the opposite sign. For the Ods/2 orbit the contribution from the triplet-even part is 
smaller than that from the triplet-odd part and they have the same sign. To see the effect 
of the tensor force on the valence proton, the energy differences between 23 F and 22 are 
shown in the table. The differences are dominated by the triplet-even part. It means that 
the contribution to the /s-splitting from the tensor force mainly comes from the triplet-even 
tensor force. 

TABLE V: /s-splitting for the proton d-orbits in 23 F with various effective interactions (see the 
text). A(Vls), A(Vt), and A(others) are the contributions to the /s-splitting from the LS force, 
the tensor force, and the other forces including the kinetic term respectively. Those are given in 



MeV. The experimental value for A(0d 3//2 - 0d 5/2 ) = 4.06MeV. 





A(0d 3/2 -0d 5/2 ) 


A(Vls) 


A(Vt) 


A (others) 


MV1 


4.2 


7.8 


-3.1 


-0.5 


MV1 without V T 


7.2 


8.3 


0.0 


-1.1 


Gogny D1S 


8.5 


9.4 


0.0 


-0.9 


M3Y-P2 


7.6 


9.2 


-0.4 


-1.2 


GT2 


8.2 


12.2 


-3.3 


-0.7 



Finally we compare the /s-splitting calculated with other effective interactions with our 
result discussed above (MV1) in Table IVl We also show the result without the tensor force 



(MV1 without Vt). The Gogny D1S force 33J does not have a tensor part and a stronger 



LS part (lTo=130MeVfm 5 ) than one we adopted above. The M3Y-P2 force [34J has a weak 



tensor part and an LS part comparable to the Gogny D1S force. The GT2 |15l. |35|| force has 
a tensor part comparable to that in the free space and a strong LS part QFo=160MeVfm 5 ). 
While the rather schematic form of the tensor force is adopted in Ref. [l5|, we replace the 
tensor part of the GT2 force with the G3RS force we used above. The sizes of the Is- 
splitting for the MV1 force without the tensor force, the Gogny force, and the M3Y-P2 force 
are large compared to the experimental value. It indicates that the relatively strong tensor 
force comparable to that in the free space is needed to reproduce the /s-splitting in 23 F. 
Although the GT2 force has a strong tensor part, it gives quite large splitting. It is due to 
the strong LS part of the GT2 force. The contribution from the LS force to the /s-splitting 
is much larger than those with other effective forces. In fact, the /s-splitting for the Op 
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orbits in 15 with the GT2 force is 8.3MeV. It is much larger than the experimental value. 
It indicates that the proper strength of the LS force, which give the reasonable /s-splitting 
in 15 is needed to reproduce the /s-splitting in 23 F. 

The tensor force also induces a 2-particle-2-hole (2p2h) correlation, which cannot be 



treated in a usual mean field calculation, 
the large attractive energy in nuclei 



36 



he 2p2h correlation by the tensor force produces 



371 ] . Recently we developed a mean field framework 



which can treat the 2p2h tensor correlation by introducing single-particle states with charge 



and parity mixing 



38 



39 



40, 



411 ] . We applied the extended mean field model to sub- 



closed-shell oxygen isotopes [41[ and found that the potential energy from the tensor force 



is comparable to that from the LS force. The im por tance o: 



the /s-splitting is indicated in other studies 



40 



42 



43 



the 2p2h tensor correlation for 



44j. It is interesting to study the 



effect of the 2p2h tensor correlation on the /s-splitting with our extended mean filed model. 
Because our calculation showed that the excess neutrons around 16 O do not contribute to 



the 2p2h tensor correlation strongly 41], the Hartree-Fock calculation seems to be sufficient 
as the first step. 



IV. SUMMARY 

We have performed the Hartree-Fock calculation with the tensor force for 15 0, 16 0, 17 0, 
22 0, and 23 F to study the effect of the tensor force on the /s-splitting. 

The tensor force does not affect the /s-splitting for the Op orbits in 15 O because 16 O is a 
LS-closed-shell nucleus. The /s-splitting is almost produced by the LS force in 15 0. 

In 22 O, the neutron 0d 5 / 2 orbit is fully occupied. It gives the finite expectation value 
for the tensor force in 22 O. In 23 F a proton is added to 22 O. The LS force works to provide 
the /s-splitting for the proton 0c?-orbits in 23 F by 7.8MeV. In contrast, the tensor force 
reduces the /s-splitting by 3.1MeV. The effect of the tensor force mainly comes from the 
occupied neutron Od^n orbit. The resulting /s-splitting of 4.2MeV close to the experimental 
data is realized by the cancelation between the effects of the LS force and the tensor force. 
The contribution from the tensor force to the /s-splitting in 23 F mainly comes from the 
triplet-even part of the tensor force. 

We have compared the results with various effective interactions with and without the 
tensor force. The effective interaction without the tensor force or with the weak tensor force 
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does not explain the experimental value for the /s-splitting for the proton Oci-orbits in 23 F. 
Our study indicates that the LS and tensor forces with reasonable strengths are needed to 
reproduce the /s-splitting in 15 and 23 F, simultaneously. 
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